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A real-time digital signal combining system for use with Ka-band feed arrays is 
proposed . The combining system attempts to compensate for signal-to-noise ratio 
(SNR) loss resulting from antenna deformations induced by gravitational and atmo- 
spheric effects . The combining weights are obtained directly from the observed sam- 
ples by using a “sliding-window” implementation of a vector maximum-MceJihood 
parameter estimator. It is shown that with averaging times of about 0.1 second, 
combining loss for a seven-element array can be limited to about 0.1 dB in a realistic 
operational environment. This result suggests that the real-time combining system 
proposed here is capable of recovering virtually ail of the signal power captured by 
the feed array, even in the presence of severe wind gusts and similar disturbances. 


I. Introduction 

There is considerable interest at the present time in 
operating the Deep Space Network at increasingly higher 
carrier frequencies for the purpose of enhancing its capa- 
bilities. To date, X-band (8.4 GHz) has been the highest 
carrier frequency employed for a deep-space mission due 
to the maturity of radio frequency (RF) components and 
technology at this frequency. Several recent studies con- 
cluded that Ka-band (32 GHz) is the proper next step 
in deep-space communications, providing 8- to 10-dB im- 
provements in downlink telemetry capability [1]. These 
potential improvements can be attributed directly to the 


higher Ka-band carrier frequencies, which yield greater an- 
tenna gains as well as reduced sensitivity to plasma ef- 
fects and increased useful bandwidth. However, there are 
some disadvantages associated with the use of higher car- 
rier frequencies, namely more stringent pointing require- 
ments, increased losses due to weather effects, and greater 
sensitivity to imperfections of the reflecting surfaces of 
the antennas [2,3]. Such imperfections become particu- 
larly troublesome on large receiving antennas, which are 
subject to significant gravitational and thermal deforma- 
tions, focussing and collimation problems, and mechani- 
cal and wind-induced vibrations. These imperfections and 
pointing errors translate directly to signal loss, resulting 
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in a corresponding degradation in telemetry performance. 
However, some of the signal can be recovered by means of 
a properly designed feed array. Our objective is to demon- 
strate a compensation technique that operates with negli- 
gible combining loss and hence recovers virtually all of the 
signal captured by the feed array. 

A conceptual design for a possible real-time antenna 
compensation system is shown in Fig. 1. Here it is assumed 
that the received signal is in the form of a temporally mod- 
ulated plane wave generated by a distant spacecraft. The 
Cassegrain receiving system consists of a large primary 
reflector, a secondary reflector, a focal-plane feed array, 
and signal processing and combining equipment designed 
to reconstruct the degraded signal. A possible feed-array 
configuration is also shown, using a maximally compact 
pattern of circular feed elements. This results in a hexag- 
onal array geometry around the central feed, with each 
succeeding “ring” adding six more elements than the pre- 
vious ring to the total; thus, one obtains 1, 7, 19, 37, etc., 
array elements as the number of rings is increased. 

The ideal primary reflector surface and its associated 
signal power distribution over the focal plane are indicated 
by solid curves in Fig. 1, while the deformed surface with 
its power distribution is shown by dashed curves. Note 
that antenna imperfections always increase the effective 
spread of the signal power in the focal plane, thus reducing 
the maximum possible signal power coupled into any single 
feed element. However, the total signal power coupled into 
a properly designed feed array may not be greatly reduced, 
although its distribution over the feeds may change with 
time. Therefore, it is reasonable to assume that degra- 
dations in telemetry performance can be ameliorated by 
the use of a matched feed array in conjuction with ap- 
propriate real-time signal processing and signal-combining 
techniques. 

Signal-processing and combining operations should be 
performed in a manner compatible with proposed DSN 
plans and equipment for future deep-space missions; the 
signal-combining system must interface properly with fu- 
ture DSN receivers. A block diagram of the proposed DSN 
“Advanced Receiver” front end is shown in Fig. 2(a). This 
receiver was designed to operate with a single feed, thus 
it has only one RF input and one complex baseband out- 
put (here represented by in-phase and quadrature com- 
ponents). The receiver front end consists of a low-noise 
amplifier, a downconverter chain that generates an inter- 
mediate frequency (IF) spectrum centered around 10 MHz, 
an automatic gain control (AGC) circuit to maintain con- 
stant average power, an 8-bit analog-to-digital (A/D) con- 
verter operating at 40 megasamples /sec, and in-phase and 


quadrature phase reference samples obtained from a digi- 
tal phase-locked loop, which downconverts the IF samples 
to baseband. The filtered outputs contain the entire mod- 
ulated spectrum shifted to baseband. 

As a first approximation, the front end of the pro- 
posed feed-array receiver is taken to be an array of receiver 
front ends (RFEs), each one associated with an element of 
the feed array. Under typical operating conditions, sev- 
eral array elements observe the signal simultaneously but 
with different amplitude and phase. Under favorable con- 
ditions, the distribution of signal power tends to change 
slowly with time as the antenna deforms due to gravita- 
tional and thermal loading. Antenna deformations and 
pointing errors introduced by wind gusts and turbulence- 
induced wavefront distortions tend to be rapidly varying 
effects, resulting in much less favorable conditions for re- 
ception. 

A block diagram of the proposed feed-array combiner 
is shown in Fig. 2(b). It consists of K separate RFEs, each 
of which generates baseband I and Q signals that serve as 
inputs to a “Signal Combining System.” Typically, K — 
1, 7, 19, and so on. This system combines the I\ complex 
inputs and generates a single complex output. The receiver 
phase-locks to the residua] carrier and provides an identical 
frequency reference to all RFEs simultaneously. Any phase 
difference between the reference signal and the received 
signal causes a measurable change in the complex output 
of the RFE. The signal combining system measures each 
complex output and uses these estimates to increase the 
effective signal level in the combined output. 


II. NASA Telemetry Format 

The telemetry format employed by deep-space ve- 
hicles has been, and presumably will continue to be, 
pulse code modulation/phase shift keying/phase modula- 
tion (PCM/PSK/PM) or PCM/PM [4]. In either case, 
the modulated carrier can be represented mathematically 
as [4] 


s T {t) = s/2I\ sin 


N 




i‘=l 


(la) 


where 


di(t) 


] PCM/PM 


* (O (0 = { (lb) 

di(t) 5\n(u> aci t) ; PCM/PSK/PM 
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6{ is the modulation index, and di(t) — il, i — 1, 2, . . . , N. 
Typically N = 1, although theoretically more than one 
data channel could be modulated simultaneously onto a 
carrier. For the case N = 1 let 8 { = 8 } in which case the 
modulated carrier becomes 

s T {i) = \/2fy[cos(6)s'm(u> c t) + s w (t)sm(6) cos(u> c t)] 

( 2 ) 

It is clear that the total received signal power is P T 
watts, the power remaining in the carrier is Prc = 
p T cos 2 (6) watts, and the power in the modulation side- 
bands is Pm = Pt ~~ Prc watts. Since for 8 = 90 deg 
there is no power in the carrier component, this type of 
modulation is referred to as “suppressed-carrier modula- 
tion,” while the term “residual-carrier modulation” is re- 
served for the case 0 < <5 < 90 deg. To date, all DSN 
deep-space probes have employed the residual-carrier for- 
mat [4]. Thus, residual-carrier modulation is assumed in 
this study, restricted to the case jV = 1 for the sake of 
simplicity. The resulting signal spectrum consists of an 
impulsive term at the carrier frequency due to the resid- 
ual carrier, and modulation sidebands centered around the 
fundamental subcarrier frequency and its harmonics. 

A graphical representation of a typical received power 
spectrum is given in Fig. 3(a), showing spectral compo- 
nents around the fundamental frequency as well as around 
the higher-order harmonics of the square- wave sub carrier. 
A typical power spectrum of the complex baseband sam- 
ples is shown in Fig. 3(b), consisting of a low-pass spectrum 
(the baseband version of the residual carrier) and the mod- 
ulation spectra around the subcarrier fundamental and its 
harmonics. Observe that the dc component can be easily 
separated from the modulated sub carrier components by 
means of a narrow-band low-pass filter, provided the fun- 
damental subcarrier frequency is significantly greater than 
the modulation bandwidth. 

The effective bandwidth of this narrow low-pass filter 
will be denoted by 5 a > whereas the effective bandwidth of 
the modulated subcarrier will be denoted by B b in sub- 
sequent analysis. In particular, it will be convenient to 
consider digital filters that perform a finite averaging op- 
eration on the input sequence, in which case these band- 
widths can be associated more precisely with the first zeros 
of the filter transfer function. 


III. System Model 

A block diagram of the signal combining system is 
shown in Fig. 4. The inputs to the combining system are 


considered to be complex samples generated by the RFEs 
associated with each feed. The sample bandwidth is as- 
sumed to be large compared to the sub carrier frequency so 
that the signal modulation is not distorted by the sampling 
operation. Since both the amplifier and background radi- 
ation effectively add noise to the signal, independent zero- 
mean complex Gaussian noise samples are added to each 
channel, with variance determined by the strength of the 
total noise process and the sample integration time. The 
noise-contaminated samples are split into two streams, one 
of these (stream A) is input to the parameter estimation 
subsystem, while the other (stream B) serves as input to 
the combining subsystem. Stream A is filtered by a bank 
of low-pass filters that removes the modulation sidebands 
in order to simplify the estimator subsystem, while stream 
B enters the combiner unperturbed. Thus, the samples 
in stream B contain both the modulation sidebands and a 
complex constant due to the residual carrier. 

It should be emphasized that the signal-combining 
system considered here consists of two separate subsys- 
tems: the parameter estimator and the channel combiner. 
The combiner structure does not depend on the form of 
the parameter estimator. Thus, various types of param- 
eter estimators could be employed to determine the re- 
quired complex weights, with varying degrees of complex- 
ity and performance. Here attention shall be restricted 
to maximum-likelihood parameter estimators, which yield 
the smallest estimation errors in the absence of a priori 
signal and noise statistics. 

The received signal in the kth channel can be modeled 
as 

s[ l \t) = V2S k [ cos(<5)sin(cj c t + 0 k ) 

-h s( l \t) sin(<5) cos (u<4 + Ok )] (3) 

where 0 k and S k = y/K are the unknown phase and am- 
plitude introduced by the antenna deformation, P k is the 
signal power in the kih channel, 8 is the modulation in- 
dex, w c is the carrier radian frequency (nominally 2n x 32 
Grad/sec) and 

$^\t) = di(*)5in(u? ac * + 0 SC ) (4) 

is the data-modulated square-wave subcarrier at a funda- 
mental frequency of lj sc / Hz. Note that the total re- 
ceived signal power Pt is the sum of the individual signal 
powers, that is, Pt = ^2 k ^ * * Both 0 k and S k are assumed 
to be slowly varying functions of time, and may be consid- 
ered constant over time intervals on the order of seconds 
to possibly minutes. The subcarrier phase 0 SC is taken to 
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be independent of the carrier phase 8 ki as it is typically 
generated by a separate oscillator on the spacecraft not 
coherent with the carrier. 1 Since d(t) and 5in(u?, c t + 0 SC ) 
both take on the values ±1 at all times, it follows that s(t) 
is also restricted to the values ±1. 

The received waveforms are downconverted to an IF 
frequency of about 10 MHz by means of a downconverter 
chain that mathematically corresponds to a single down- 
conversion operation. It is assumed here that downcon- 
version is performed with a “perfect” frequency reference 
that does not introduce any additional phase errors into 
the resulting IF signals. These signals can be represented 
as 


Here it is assumed that the sampling time T 0 is small 
enough to avoid aliasing and that there are sufficient quan- 
tization levels to allow the representation of the quantized- 
amplitude samples by continuous-amplitude samples. 

The final down conversion to baseband is accomplished 
using digital samples. Both in-phase and quadrature sam- 
ples are generated simultaneously using samples of a per- 
fect IF frequency reference, as indicated in Fig. 2(a). This 
operation is represented mathematically as 

sSk(i) - [2«/fc(i)sin(w/fT 0 )] LP (8a) 

cSk(i ) = [2s/jfc(*’) cos(u>/z7o)] Li > (8b) 


4V (0 = [4 1} (0>/2 cos ((w e - w/)o] Lp 

= Sk jcos(<5) sin(w/t -j- $k ) 

+ sin(5) cos(ujjt -{- Ok)] (5) 

where “ LP ” refers to a low-pass filtering operation. 

The IF waveforms are sampled by an A/D converter 
which generates 8-bit samples every T 0 seconds, or in the 
case of the Advanced Receiver, at a rate of 40 million sam- 
ples per second. It is convenient to view these samples 
as the coefficients of an orthonormal expansion with basis 
functions the well-known “sampling functions” 


The sample products are filtered by digital low-pass filters 
that pass only the modulated subcarrier frequencies. The 
resulting baseband in-phase and quadrature samples in the 
kth channel at time iTq become 

8 Sk(i) = Sk [cos(5) cos(0*) - $(i) sin(6) sin(0*)] (9a) 

and 

c$*(z) = S k [cos(£)sin(0*) + s(i) sin(<5) cos(0*)] (9b) 

where 

s(t) = d(*To)5in(w fe *To + 0, c ) (9c) 




sin [{*fT 0 )(t - iT 0 )] 
(7t/To)(<-iTo) 


( 6 ) 


which form a complete orthonormal set in the space of 
band-limited functions over the real line. In terms of these 
basis functions the expansion becomes 


CO 

4V(0 = X! (7a) 

< = — 00 


These signals contain the unknown amplitude and phase 
that affect both the residual-carrier and the moduiated- 
sub carrier signals. 

Making use of the identities 

Re{e J *[a + jb]} = acos(0) - 6sin(0) (10a) 

Im{e ; *[a -f jb]} = a sin(0) + bcos(6) (10b) 


with coefficients that represent samples of the waveform at 
the sampling instants “*Tq” : 

= f 4V(0<M0 < * 

J — OO 

= St [cos(<5) sin(w/i7o + Ok) 

+ s (1) (iT 0 ) sin(5) cos(w/iT 0 + 0 k )] (7b) 


1 P. Kinman, Telecommunications Systems Section, personal 
communication. 


and letting “ denote a complex quantity, one can also write 
,«*(») = Re{s*(i')} (11a) 

c Sk(i) = Im{st(i)} (lib) 

where 

Sk(i) = Ske j9k [cos(£) + js(i) sin(5)] 

= V* [cos(<5) + js(i) sin(<5)] (12) 
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is the “complex envelope” of the baseband signal samples. 
The complex coefficient V k contains both the amplitude 
and phase of the baseband sequences. 

Next the representation of the additive noise is considered. 
Bandpass noise in the kth channel can be modeled as 

n k {t) = \/2[n c t(0 cos {u e t) + n, k (t) sin(w c <)] (13) 

where n ck (t) and n, k (t) are independent zero-mean Gaus- 
sian baseband processes, each with spectral level No k /2. 

Downconverting to IF yields 

n n(t)= [ntCOV^cosftwc -ui)t}] LP 

= n ct (<) cos(w/t) + n,jt(0sin(w/<) (14) 

while the IF samples are again the coefficients of the or- 
thonormal expansion 

OO 

njk(t) = (15a) 

i= — oo 

/ co 

-OO 

= n ck (i) cos (w/iTo) + n sk (i) sin(w/iT 0 ) (15b) 

Assuming that the noise processes are broadband 
compared to the sampling bandwidth, the correlation func- 
tion of the noise samples becomes 

£’[n c t(i)n c jfc(^)] = £[n»*(O n »*M] 

= | ( 16 ) 

In-phase and quadrature noise samples are generated 
along with the signal samples, with the result that 

,n k (i) = [2n/fc(i)sin(w/»To)] iP = n, k {i) (17a) 

and 

c n k (i) = [2n Ik (i) cos(u)/tT 0 )j Lp - n ck (i) (17b) 

A complex representation for the noise samples is obtained 
by letting 

,n t (t) = Re{n t («)} ( 18 a) 

£«*(«’) = Im{ht(i)} ( 18 b) 


where 

n*(i) = n lk (i) + jn ck (i) (18c) 

is a sequence of complex noise samples. Note that the 
variance of the complex noise samples is 2crjj!, because the 
real and imaginary components are independent zero-mean 
processes, each with variance cr^. The use of complex no- 
tation will often allow considerable simplification in the 
derivations that follow. 

Next, the problem is generalized somewhat by sam- 
pling of the primitive samples in streams A and B at dif- 
ferent rates. The final samples are formed by averaging 
consecutive primitive samples, and hence represent an ef- 
fective low-pass filtering operation. Since the samples in 
stream B must reflect the temporal variations in the mod- 
ulated subcarrier accurately, typically only a few primi- 
tive samples can be averaged here. In stream A, however, 
enough samples must be averaged to ensure that the mod- 
ulated subcarrier terms are filtered out, leaving only the 
residual carrier terms. Here one must be careful to make 
sure that the underlying processes of interest are not fil- 
tered out. Therefore, in stream A strict lower and upper 
bounds must be observed on the number of primitive sam- 
ples averaged, while in stream B only the upper bound is 
relevant. 

Suppose that M* primitive samples are averaged in 
every channel of stream A, and Mb in stream B, This 
means that the time delay between samples is Ta = MaTq 
in stream A and Tb — AffjTo in stream B. The situation is 
illustrated in Fig. 5, which shows the timing relationship 
between the samples in the two streams. 

Applying the averaging operation of stream B to the 
primitive signal samples again yields Eq. (9), but evaluated 
at times is — iMsTo* The averaging operation in stream 
A effectively removes the sub carrier terms from the prim- 
itive samples, leaving the residual-carrier samples defined 
at times iA — iA/^To: 

,x k (i A ) = [» s k(i)] LP = S k cos(5) cos(0*) (19a) 

c x k (i A ) = [<:**(*)] ip = S k cos(6)sin(0*) (19b) 

Equation (19) can be expressed in complex form as 

it(»x) = ,*k(i A ) + 

= S k e? 9k cos(8) 

= Vjt cos(6) = X k = X k R -f jX k T (20) 
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Note that samples in both streams are obtained from 
the same sequence of primitive samples in the correspond- 
ing channels, hence the resulting sample variances in the 
corresponding channels are always related by the ratio 
of the number of primitive samples averaged in the two 
streams. Let the noise samples in stream A be denoted 
by m k (i A ) 9 and in stream B by n k (i B ). The averaging 
operation reduces the noise-sample variances by factors 
of Ma and M Bl respectively, so that ff Ak = a I/M a and 
a Bk = a \IM B . The noise variances in the two streams are 
related by <j\ k = r}<r 7 Ak , where rj = M A /M B is the ratio of 
the effective band widths of the two streams. 

The parameter-estimation problem is greatly simpli- 
fied if estimation is based only on the residual-carrier se- 
quence as given by Eq. (20), because interference from 
the data modulation is not present. Since only a small 
fraction of the total signal power is in the residual car- 
rier, one might expect substantial performance improve- 
ments by making use of the modulation sidebands as well. 
However, the more complicated estimation algorithms that 
make use of the full-spectrum sequence of Eq. (12) will not 
be addressed here. It is emphasized that while the averag- 
ing operation in stream B is optional, it is mandatory in 
stream A whenever modulation sidebands are present. 


Although in general p z is a function of time, here it is 
assumed that all relevant processes are stationary, so that 
the time index can be ignored. Substituting Eq. (21b) for 
?*(**) yields 


Pz = 


E*=i [cos(6) + j sin(<5)s(j' B )] 


2Zi = l 1^*1 2<Tg k 


(24) 


Making use of the Schwarz inequality, 
V k 


K ( V > 

X] (V2u)i 0-flit) — [cos($) + j sin(<5)s(i B )] 


K K Vt 

<^2\w k \ 2 2a 2 Bk J2^ 

k = l ■ • 


(25) 


k = 1 Bk 


Dividing both sides of Eq. (25) by Ylk \™k\ 2 2v Bk yields 
the following upper bound on p 2 : 


K 


Vtl 


i = l 2<T Bk 


A 

= P 


(26) 


IV. The Combining Algorithm 

The mathematical model for the combining system 
depicted in Fig. 4 can be summarized as follows: the noise- 
corrupted complex baseband samples in streams A and B 
are modeled as 


The choice of weights that achieves the value of the right- 
hand side of Eq. (26) is 

"* = it (2?a) 


fijfc ( M ) = £*(*’>0 + ^ At ( m ) 


(21a) as 'can be seen by direct substitution 


h(iB) = s k (i B ) + n k (i B ) (21b) 

Complex combining weights are determined by the esti- 
mator subsystem operating on stream A, and applied to 
the samples f*^) in stream B to produce the combined 
output sequence z(i B ): 


K 

z{ib) = y^r k (i B )w k 

k=l 


( 22 ) 


The object of the combining algorithm is to maximize the 
signed- to-noise ratio of the combined output samples, de- 
fined as 


a |g[f(»B)]| 2 _ 

- var [*(**)] " var[Ef =1 ^f*( lB ) 


(23) 



K 


= E 


I vS 


k = 1 2<T Bk 


-p 


(27b) 


Therefore, if the signal and noise parameters were known, 
then the combining operation defined by Eqs. (22) and 
(27a) would achieve the sample SNR p defined in Eq. (26). 
Note that with the optimum choice of weights, the mean 
value of the combined signal and the variance of the com- 
bined noise are both exactly equal to p. If all noise vari- 
ances are equal, a Bk ~ a B , then we have the interpretation 
P — Pt where Pt = J2k 15 total signal power 
collected by the array. 
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Since the required signal and noise parameters are 
generally not known, this upper bound on the sample SNR 
is not achieved in practice. When dealing with unknown 
parameters, a plausible approach is to use the best avail- 
able estimates of the underlying parameters to estimate 
the combining weights. This is the approach adopted in 
this article; the estimated weights are obtained by applying 
the optimum weight formula (Eq. 27a) to the maximum- 
likelihood estimates of the underlying parameters {Vk} and 

V. The Maximum-Likelihood Estimator 

The combining algorithm requires an estimate of the 
complex signal coefficient Vk and of the noise variance a Dk • 
These quantities can be estimated simultaneously from the 
observed samples, either in stream A or in stream B. Here 
attention is restricted to the residual-carrier sample stream 
(stream A), because of the relative simplicity of the prob- 
lem when the data- modulated subcarrier is not present. 

In the absence of a-priori statistics about the values 
of the unknown parameters X* and the maximum- 

likelihood approach yields the best estimates of the em- 
bedded parameters, on the average. These estimates are 
obtained from the likelihood function, which is the con- 
ditional joint probability density of the observed sam- 
ples, conditioned on the parameters Xk and o\ k * The 
maximum-likelihood estimates are those values Xk and 
cr Bk that simultaneously maximize the likelihood function, 
or a monotonically increasing functional of it. 

The log-likelihood function is obtained from the con- 
ditional probability density of the observables, conditioned 
on the parameters of interest. In stream A only the resid- 
ual carrier is present, hence the observed complex samples 
in the Jbth channel at time *i A T A n can be represented as 

ujb(iyt) = Xk + rh k (i A ) (28) 


where X k = V* cos(5), and rh k (i A ) is a zero-mean com- 
plex Gaussian random variable with variance 2cr^ fc . It is 
assumed that all noise samples are independent random 
variables. If L consecutive samples are observed in each 
of K channels, the joint probability density of the entire 
K x L complex sir ray becomes 


K L 


P (u|x,<ri)=n n ( 2 *du) 


2^-ig-l n k {i A )-X k \ 2 !2a\ k 


k=li A =l 


(29) 


where u is a K x L-dimensional complex matrix whose el- 
ements encompass all possible values of the indices i A 
and “jfc ,” while X and <t\ are X-dimensional complex vec- 
tors. Denoting the real and imaginary components of X k 
by X k R and A*/, respectively, the natural logarithm of 
Eq. (29) becomes 

A = ln[p(«|X,*i)] 

K 

= -L'E,M**°Ak) 

k=l 

t^A Z<T Ak iA=1 

+ («*/(;*) -x*/) 2 ]^ ( 3 °) 

This is called the log-likelihood function. The maximum- 
likelihood estimates of X k R, X\ t/ , and (r\ k are those values 
X kR) X kI} and b\ k that simultaneously maximize the log- 
likelihood function A. Differentiating A with respect to 
each XkR and X*/, setting the result equal to zero, and 
solving yields the maximum-likelihood estimates 

XkR = £ ^2 

X k i = y (3it>) 

L 1=1 

and 

X k — XkR + jXki k = 1 1 2, . . . , K (3 1 c) 

The estimate of the noise variance is obtained by writing 
the log-likelihood function, maximized over X, as 

T — < 32 ‘> 

fc=l h = l\ Ak / 

G k 4 £ [(«w(m) - + («*/(*'>») - *«)’] 

' A = 1 (32b) 

Differentiating with respect to °Ak , setting the result equal 
to zero, and solving yields 

•i* = 5Z°‘ <33) 
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Note that the maximum-likelihood signal estimates 
XkR and X\ b/ are unbiased estimates of the mean values 
of the real sequences UkR{i A ) and Uki(i A ), also known in 
statistics as the “sample mean.” However, the maximum- 
likelihood estimates of the noise variances are biased, and 
hence do not correspond exactly to the unbiased “sam- 
ple variance” used in statistics. Here we use the unbiased 
variance estimate 


-2 

a Ak t ub 


1 


2(L-1) 


G k 


(34) 


<*Ak,ub(*A) — &Ak,ub(iA ~ 1) + — jy — 1)| 2 

- («t(M - L - 1)| 2 

- L [|^*(m)| 2 - |i*(u - 1)| 2 ) (36b) 


Note that the corrections to the previous estimates can be 
obtained from a digital “first-in, first-out” stack, minimiz- 
ing the number of computations required for each update. 


For large L, the difference between the two estimates be- 
comes negligibly small. The unbiased estimates shall be 
used in this analysis because each of these estimates can 
be modeled conveniently as the sum of the true variance 
plus a zero-mean error term. 

The maximum-likelihood estimators defined above 
provide the best possible estimates of the desired param- 
eters if a-priori statistics are not available, as long as the 
signal parameters and noise statistics do not change with 
time. These results can be extended to the slowly varying 
case by selecting the total observation time to be small 
compared to the characteristic time scale of the variations. 
Thus, over a short time interval the above solutions apply. 
However, the observation interval must now be continu- 
ously shifted in time to follow the parameter variations. 
This extension leads to the following “sliding-window” es- 
timator structures: 


*A~ 1 


Xk(i A ) = t £ (35a) 

1=ia-L 

1 ^ 1 m 2 
.uiO'iO ~ 2(17- 1) ^ (35b) 

? t — i a~ L 


Thus, the current estimates at time “i A T A ” are based on 
the previous L samples, up to and including the sample at 
time 

For large values of L and short sample durations, the 
computational burden required by these estimators may 
become quite severe. This difficulty can be ameliorated 
by observing that a recursive implementation is possible, 
where the current estimate is obtained in terms of the pre- 
vious estimate, the current observable, and the “oldest” 
observable in the previous sum. The recursive forms are 

X k (i A ) =X k (i A - 1) + i [u k (i A - 1) - u k (i A - L - 1)] 

(36 a) 


It is well known that the maximum- likelihood estimate 
of the ratio of two parameters is equivalent to the ratio 
of the maximum-likelihood estimates [5]. It follows that 
the maximum- likelihood estimate of the £th weight can be 
obtained from the maximum-likelihood estimates of the 
complex signal estimates and the noise variances as 


^Jt(M) = 


V k (i A ) _ X k (i A ) 

2w 2 Ak (i A ) 2cos(6)i]a 2 Ak (i A ) 


(37) 


The maximum-likelihood estimate of the Jbth weight there- 
fore depends on the signal parameters in each channel of 
stream A and on the corresponding noise variances. 

An interesting special case occurs if the noise variance 
in each channel is the same. This situation could arise if 
the amplifier characteristics in the various channels were 
exceptionally well matched and all feeds observed the same 
background power levels, so that the noise temperatures in 
the various channels were essentially the same. In this case, 
the common variance does not have to be estimated, since 
it becomes a scaling factor applied to all of the weights. 
Because each signal and noise term is equally scaled, the 
SNR of the combined noisy signal is not affected if this 
common scaling factor is ignored. 


VI. Combiner Performance Analysis 

In this section, the performance of the signal combin- 
ing system will be evaluated. Constant signal components 
will be assumed throughout, and it shall also be assumed 
that L is large, so that the biased and unbiased variance 
estimates can be considered equal. However, the unbi- 
ased variance estimates will be used in the analysis. Per- 
formance will be evaluated in the “steady-state,” that is, 
assuming that L or more samples have already been ob- 
served. 

Assuming the “sliding-window” configuration of 
Eq. (35), it is seen that the estimates at time i A T A are 
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based on the previous L samples, up to and including the 
sample at time (i A - 1)7Y Since the estimates are un- 
biased, each can be expressed as the sum of the true pa- 
rameter plus a zero-mean random variable representing the 
estimation error. This model yields 

X k (i A ) = x k + (Z k (i A ) - Xk) = x k + € Xh (iA) (38a) 
&Ak(iA) = 0Ak + {&Ak(*A) ~ a Ak ) = °Ak+ f <>%,,(**) (38b) 


3 <j\ k . Since for a Gaussian population the sample mean 
and sample variance are independent random variables, it 
follows that ?x(0 an d to 7 are independent as well. Since 
the i A th estimate is basecf on L previous samples up to and 
including sample number (i A — 1), and since each sample 
is independent of all others, it also follows that the i A th 
estimate is independent of the i A th sample. 

Making use of the model defined in Eq. (38), the 
maximum-likelihood estimate of the fcth weight becomes 


where 

^x fc (u) = ^ 53 (39a) 


W k {tA) = 


XI 

2 cos(<5)cr 


2 

Bk 


(l+gX fc (»A)) 

( 1 + «*i*(*^)). 


(41a) 






Y |a*(o 

(39b) 


where 


= ?xJXk 


(41b) 


Since the estimates are unbiased, it follows that the error 
terms are zero mean. The variance of the estimation errors 
can be obtained using well-known expressions for the vari- 
ance of the sample mean and the variance of the sample 
variance, as follows: 


var[e Xk (u)] = 


E\rh k {i A )\ 2 

L 


^f\k _ ^ cr lt 

L r]L 


(40a) 


and 

e °l* = € *\,J°Bk ( 41C ) 

If the standard deviation of each estimatiori error is small 
compared to the estimate, then Eq. (41a) can be approxi- 
mated as 


var 


= 5 var { (T“T) L \ UkR{e) ~ } 


-s(' 


1-3 4 

2 L ~ ~ 1° M 


’Ak 


L- 1 


(40b) 


™ k{iA) = 2cosW47 [1 + 

x [i + 4| „(*'*)-•••] 


x; 


2 cos (6)<r# k 


[l + “ e <^» (**))] 


XI 


2 cos (S)(Tg k 


[l + e^(u)] 


while 

var kM = rrf = r^r (40c) 

In deriving Eq. (40b) recall the fact that the fourth 
central moment of a Gaussian random variable is = 


= + 1 < 42 > 

za Bk 

Therefore, the weight estimates may also be expressed 
approximately as the sum of the true weight plus an error 
term, provided the errors in the component estimates are 
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sufficiently small. As before, the error term is zero mean, 
with variance 

var K*(u)] = var [^(m)] + var [^^(u)] 

1 




+ 


< 2 (l - 1) 


x k L 




n v k 


+ 


L cos 2 ((5) 


(L-l) 


(43) 


var [s c (i'u)] = 


»/Lcos 2 (6) ( L 


1 K 


! ) fc 4 ^ k 


(47) 


Note the dependence of the signal variance on the sum of 
squares of the channel SNRs in Eq. (47). This term can 
be bounded as follows: 


2 K 

-<y 

K - A& 


_ 4 

Vk 


<p l 


*=1 


Bk 


(48) 


Note that the current weight estimate is independent of 
the current sample, since it depends only on the previous 
L samples. Using the “small-error” approximation for the 
weighting factors, the combiner output sequence defined in 
Eq. (22) can be expressed as 

z(i B ) ^ s c (i B ) + n c (i B ) (44) 

where 


The left side follows from the Schwarz inequality: equal- 
ity holds when all the channel SNRs are equal. The right 
side is trivial since one obtains this inequality by throwing 
away all of the cross terms in the square of the series for p . 
Equality holds when all the signal power is in a single chan- 
nel. Useful bounds on p follow if the noise spectral levels 
in the various channels are bounded by known minimum 
and maximum values No min and No ma r : 


K v k \ 

Sc(*b) = [cos((5) + j sin(<5)s(i B )] E [l + e^fo)] 

*=i la B k 


(tt — ) ToMb =Pi<P<Pu = ( Jp—) To Mb 

\ Nq max) \ No min J 

(49) 


and 


(45a) With the help of the above bounds, the variance of the 
signal term can be bounded as 


Mia) = E Vk 'y (tB) [1 + CM (45b) 

k=l Bk 

In these equations, the index i a refers to the last estimate 
in stream A that occurred prior to the current sample i B 
in stream B. 

The combiner output in Eq. (44) can be viewed as 
having two components: a signal term with random mag- 
nitude that takes on a new value each time the weight index 
changes, plus an equivalent additive-noise term that takes 
on a different value for each new sample of the combined 
signal. The magnitude of the mean value of the signal term 
is p, the sum of the channel SNRs that would be obtained 
if the weights were known without error: 


& , p! 

t]Lzo& 2 (6) K{L— 1) 


< var [ s c (i B )] 


i+trbj' 50 ' 

In addition to the known system parameters, these bounds 
involve only the total received signal power and the bounds 
on the noise spectral levels. 

The equivalent noise term defined in Eq. (45b) consists 
of sums of random variables obtained from all K channels. 
Since the random processes in the various channels are 
independent, and since the current noise sample is inde- 
pendent of the current weight error, one can express the 
variance of the total complex noise as 


K 


p=y 


v k 

2a% k 


(46) 


The variance of the signal term can be expressed in 
terms of p as 


var [n c (i B )] = p 



1 K 

L-l + f/Lpcos 2 ^) 


(51) 


Note that the variance of the combined noise is always 
at least as great as the variance of the ideally weighted 
noise terms, p, approaching that limit as L grows without 
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bound. The excess noise terms can be attributed directly 
to uncertainty in the weight estimates. 

Perhaps the most direct measure of combiner per- 
formance is the extent to which the sample SNR of the 
combined signal approaches its maximum value, p, which 
would be obtained if the complex combining weights could 
be determined without error. A lower value here implies 
degraded performance. Analogous to the previous defini- 
tion in Eq. (23), the sample SNR of the combined signal 
when the maximum-likelihood weight estimates are used 
is defined as 


_ 1 £[*-(»b)]I 2 

Pml var [5 (*b)] 


(52) 


The time index on pml is ignored since all relevant 
processes in this model are assumed to be stationary. The 
magnitude of the expected value of the combiner output 
is again p, while the variance is the sum of the signal and 
noise variances defined in Eqs. (47) and (51): 


var [z(»b)1 = P 


1 + 


1 1 

L - 1 t)Lcos 2 (6) 




+ 


1 

p(L - 1) 


K 


£ 


Vi 

4 <T B* 


(53) 


7» = 


1 K 


m 


t)L cos 2 (6) p(L - 1) j“ 4 ^Bifc 

1 a 

~ ^ ^ L — 1 rjLpcos 2 (6) 


(55b) 

(55c) 


Since these loss factors involve both sums of the individual 
sample SNRs and sums of their squares, it is useful to 
invoke the bounds of Eqs. (48) and (49) here. Thus, the 
component loss factors are bounded as 


A 

a* = 


pi 


r)L cos 2 (<$) ^ K{L- 1) 


<7. 


1 Pu A 

~ i]Lcos 2 (6) (L — 1) 


(56a) 


*7n =1 + 

< 1 + 


K 


+ 


1 


r}Lp u cos 2 (<5) ( L — 1) 

K 


< 7n 


rjLpt cos 2 ((5) (L- 1) 


+ / r 1 — = u7 n (56b) 


while the total loss is loosely bounded by sums of the upper 
and lower bounds on its components: 


n» + tin = tl ^ 7 ^ uT — u7i + njn (56c) 


Substitution into Eq. (52) yields 
Pml — 


P 

+ ITT + nLcL'Wj + 7 ) + p(i-l) 4^7 

(54) 

Note that pml is always less than p , approaching that 
maximum value as L approaches infinity. This is reason- 
able, since the maximum-likelihood weight estimates ap- 
proach the true weights as the number of observed samples 
increases. A direct measure of the loss in sample SNR due 
to the imperfect weight estimates is the quantity 


Note that tighter bounds on the total combining loss 
exist, but these involve complicated formulas in different 
ranges of the parameters, and will not be considered here. 

The equal noise-variance case can be treated as the 
general case considered above, but with considerable sim- 
plifications throughout. Although the noise variance is not 
estimated in this case, a common noise variance denoted 
by <t 2 b shall be included in these equations for consistant 
results. Now the combining weights are given by 


U'i(jyl) = 


£(U) 

2 cos(6)<Tg 


(57) 


where the denominator contains but not its estimate. 
The weight estimate can now be expressed as 


7 — p! Pml 

- J, + In 


v: 


(55a) 


“>*(**) = 2^2 l 1 + 


(58) 


with components due to signal, 7 ,, and noise, j n) defined 
as 


and it is seen that there is no need to make a small error 
approximation here. Therefore, these results will be valid 
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in general, not only where the small-error approximation 
holds. The variance of the combined signal and of the 
combined noise now becomes 


var [s c (*b)] 



( 1 
\r}L cos 2 (6) 


) 


= p 
r)L cos 2 (S) 


(59a) 


var [n c (? 5 )] - p 



K 

rfLp cos 2 (<5)_ 


(59b) 


Hz. Assume a seven-element array, K = 7. As before, let 
To — 2.5 x 10 -8 sec, corresponding to a first-zero frequency 
of /o = 1/T 0 = 4 x 10 7 Hz. Since the modulated subcar- 
rier must be filtered out of stream A, the first zero of this 
effective filter should be at some frequency much smaller 
than the subcarrier fundamental, but large enough to pass 
the slowly-varying parameters generated by the antenna 
deformations. This requirement yields f A < 10 4 Hz or so, 
implying that M A = fo/f A > 4 x 10 3 . Observe that near 
the above parameter values, the “small-error” approxima- 
tion upon which the general formulas are based remains 
valid, provided the inequality L 1 is satisfied. 


while the variance of the total combined samples is their 
sum, as before: 


var [z(*b)] = var [s c (*b)J 4- var [n c (*s)] 

(>+*) 


- P 


1 + 


i }L cos 2 (<5) 


(59c) 


The combining loss for the equal- variance case also follows 
directly 


tl = tin 4 els 


(60a) 


with components 


tin — 1 4 * 


tls = 


I< 

r]Lpcos 2 (S) 

1 

r\L cos 2 (<5) 


(60b) 

(60c) 


The behavior of the total combining loss and of its com- 
ponents will be examined in the next section. 


VII. Numerical Results 

The following example serves to illustrate the behav- 
ior of the combining loss in a realistic setting. Suppose 
the 70- m antenna at DSS 14 is receiving signals from a 
deep-space vehicle at Ka-band, and further suppose that 
the received power levels and noise temperatures are sim- 
ilar to those expected from the Galileo spacecraft during 
its encounter with Jupiter, that is, Pt/Nq = 55 dB-Hz 
= 3.16 x 10 5 Hz. Also assume that a standard NASA 
modulation format has been employed with 5 = 80 deg, 
sub carrier frequency of fic = 2 x 10 5 Hz, and assume that 
we wish to recover the fifth sub carrier harmonic (the sec- 
ond and fourth harmonics are zero for a square wave) so 
that it is required that fs > 10 6 Hz; let fs = 2 x 10 6 


The bounds on the combining loss 7 are shown in 
Fig. 6, for nominal Pr/N 0 of 45, 55, and 65 dB-Hz. Note 
that the combining-loss bounds, expressed in dB, are plot- 
ted on a logarithmic scale to show their relative behavior 
in regions where each combining-loss bound is very close 
to zero dB. The maximum and minimum values of Pt/Nq 
are taken to be ± 2 dB above and below the nominal val- 
ues presented in the graphs. Since the total averaging time 
used for each weight estimate T can be related to the num- 
ber of observed samples L as T = M A T 0 L t it is clear that 
the combining-loss bounds can be viewed either as func- 
tions of L or as functions of T. Therefore, the loss bounds 
are plotted in terms of both of these fundamental vari- 
ables in Fig. 6. Because of the small-error approximation 
the bounds may not be accurate in the region L < 10 (or 
equivalently, T < 10“ 3 sec), however, this region is of no 
interest in the present application, where the timescale of 
interest is typically on the order of a second. Note that 
a total combining loss of less than 0.1 dB can be guaran- 
teed with observation times of less than 0.1 second, when 
operating at 55 dB-Hz. This result shows that during a 
typical encounter, this combining system can operate with 
negligible combining loss, even in the presence of severe 
wind gusts which can induce rapidly varying mechanical 
distortions into the antenna structure. 

The components of the combining- loss bounds due to 
signal and noise are shown in Figs. 7(a) and 7(b), again 
expressed in dB, at a nominal Pt/Nq of 65 dB-Hz, with 
all the other system parameters the same as above. Here 
the total combining-loss bounds and their components (in 
dB) are displayed on a linear scale in order to show the 
signal-loss bounds on the same graph. The loss compo- 
nents displayed in Fig. 7 indicate that for large L the great- 
est contribution to the total loss comes from the equivalent 
additive noise. The losses due to the random variations in 
the signal tend to be relatively insignificant for large L. 
This behavior is even more pronounced at lower Pt/Nq , 
where the individual components of the combining loss are 
even more difficult to separate. However, the reader is 
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cautioned that the total combining loss defined here may 
not be an adequate indicator of performance in all applica- 
tions. The relative contribution of these loss components 
to the performance degradation of specific systems may 
difTer markedly and should be carefully examined for each 
application. 

The combining loss for the equal- variance case can be 
determined exactly, without resorting to small-error ap- 
proximations. These results are shown in Fig. 8, which con- 
firms the previous conclusion that total combining losses 
can be limited to 0.1 dB with averaging times of less than 
0.1 seconds when operating at 55 dB-Hz. The behavior of 
the loss components shown in Fig. 9 again indicates that 
the dominant contribution comes from the additive noise, 
so that losses due to random signal variations are relatively 
insignificant in this region. 


VIII. Summary and Conclusions 

A real-time signal combining system for use with ar- 
ray feeds has been proposed and evaluated. The combining 
system can be used to compensate for losses resulting from 
time-varying antenna distortions due to gravitational and 
atmospheric effects, provided the characteristic timescale 
of the variations is not too rapid — perhaps on the order of 
a second or so. It is ideally suited for recovering losses re- 
sulting from main-reflector gravity- and wind-induced de- 
formations of large DSN antennas operating at Ka-band 
or higher frequencies. 

Real-time compensation is achieved by forming a 
weighted combination of the array-feed signals observed in 
the presence of additive noise with combining weights con- 
tinuously adjusted to maximize the signal- to-noise ratio of 
the combined signal. In the digital baseband model consid- 
ered, the weights were complex numbers obtained from a 
“maximum-likelihood” weight estimator that obtained the 
required signal and noise parameters from observations of 
the array-feed outputs. A “sliding-window” implementa- 


tion was proposed that effectively follows the variations 
in the signal and noise parameters, achieving automatic 
compensation in the combined output. 

A mathematical model of the real-time compensation 
system was developed. Performance measures were defined 
and system performance evaluated under a “small-error” 
approximation, which is valid in the regions of interest. 
The model assumed constant signal and noise parameters, 
but remains valid for time-varying parameters, provided 
the variations are slow compared to the time scale of the 
observations. It was found that the combining operation 
produced two distinct effects: it introduced random varia- 
tions in the combined signal and it increased the effective 
power of the additive noise. Both of these effects could 
be attributed to the uncertainty in the weight estimates 
due to the presence of additive noise in the system. It was 
shown that under receiving conditions similar to those ex- 
pected during the Galileo encounter, the total combining 
loss could be limited to 0.1 dB, with averaging times of 
less than a hundred milliseconds for a seven-element array. 
This means that virtually all of the signal power collected 
by the feed array can be recovered, even in the presence of 
wind gusts that may produce changes in the signal param- 
eters on the timescale of a second. For larger arrays, the 
combining loss tends to increase. 

Although two different types of degradations have 
been defined and quantified, their effect on telemetry sym- 
bol detection, carrier tracking loop behavior, mechanical 
pointing-system behavior, etc., has not been completely 
determined here because the interaction of these compo- 
nents tends to be applications-dependent. These issues 
should be addressed in the future. Also, mathematical 
models should be developed to characterize the problem in 
a dynamic environment and the effect of imperfect phase 
reference on combiner performance should be addressed. 
In addition, the possibility of using the entire modulated 
spectrum for determining the combining weights should be 
addressed, since this approach will improve combiner per- 
formance in the low-SNR or highly dynamic environments. 
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SAMPLES 


Fig. 1. Real-time antenna-compensation system conceptual design. 
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Fig. 3. NASA telemetry spectra: (a) signal spectrum after IF downconverslon and (b) spectrum 
of baseband samples. 
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Fig. 4. Signal-combining system (baseband model) block diagram. 













